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Abstract — We propose a branch flow model for the anal- 
ysis and optimization of mesh as well as radial networks. 
The model leads to a new approach to solving optimal 
power flow (OPF) that consists of two relaxation steps. 
The first step eliminates the voltage and current angles 
and the second step approximates the resulting problem 
by a conic program that can be solved efficiently. For radial 
networks, we prove that both relaxation steps are always 
exact, provided there are no upper bounds on loads. For 
mesh networks, the conic relaxation is always exact but the 
angle relaxation may not be exact, and we provide a simple 
way to determine if a relaxed solution is globally optimal. 
We propose convexification of mesh networks using phase 
shifters so that OPF for the convexified network can always 
be solved efficiently for an optimal solution. We prove 
that convexification requires phase shifters only outside a 
spanning tree of the network and their placement depends 
only on network topology, not on power flows, generation, 
loads, or operating constraints. Part I introduces our 
branch flow model, explains the two relaxation steps, and 
proves the conditions for exact relaxation. Part II describes 
convexification of mesh networks, and presents simulation 
results. 

I. Introduction 

A. Motivation 

The bus injection model is the standard model for 
power flow analysis and optimization. It focuses on nodal 
variables such as voltages, current and power injections 
and does not directly deal with power flows on individual 
branches. Instead of nodal variables, the branch flow 
model focuses on currents and powers on the branches. 
It has been used mainly for modeling distribution cir- 
cuits which tend to be radial, but has received far less 
attention. In this paper, we advocate the use of branch 
flow model for both radial and mesh networks, and 
demonstrate how it can be used for optimizing the design 
and operation of power systems. 

One of the motivations for our work is the optimal 
power flow (OPF) problem. OPF seeks to optimize a 

A preliminary and abridged version will appear in [1]. 



certain objective function, such as power loss, gener- 
ation cost and/or user utilities, subject to Kirchhoff's 
laws, power balance as well as capacity, stability and 
contingency constraints on the voltages and power flows. 
There has been a great deal of research on OPF since 
Carpentier's first formulation in 1962 [2]; surveys can 
be found in, e.g., [3]-[7]. OPF is generally nonconvex 
and NP-hard, and a large number of optimization algo- 
rithms and relaxations have been proposed. A popular 
approximation is the DC power flow problem, which 
is a linearization and therefore easy to solve, e.g. [8]- 
[11]. An important observation was made in [12], [13] 
that the full AC OPF can be formulated as a quadrat- 
ically constrained quadratic program and therefore can 
be approximated by a semidefinite program. While this 
approach is illustrated in [12], [13] on several IEEE 
test systems using an interior-point method, whether or 
when the semidefinite relaxation will turn out to be 
exact is not studied. Instead of solving the OPF problem 
directly, [14] proposes to solve its convex Lagrangian 
dual problem and gives a sufficient condition that must 
be satisfied by a dual solution for the duality gap to be 
zero and for an optimal OPF solution to be recoverable. 
Importantly, this provides a way to determine for sure if a 
candidate solution is globally optimal for the nonconvex 
OPF problem. This result is extended in [15] to include 
other variables and constraints and in [16] to exploit 
the sparsity of power networks and phase shifters for 
the convexification of OPF. In [17], [18], it is proved 
that the sufficient condition of [14] always holds for 
a radial (tree) network, provided the bounds on the 
power flows satisfy a simple pattern. See also [19] for a 
generalization. These results confirm that radial networks 
are computationally much simpler. This is important as 
most distribution systems are radial. 

The limitation of semidefinite relaxation for OPF is 
studied in [20] using mesh networks with 3, 5, and 
7 buses. They show that as a line-flow constraint is 
tightened, the sufficient condition in [14] fails to hold 
for these examples and the duality gap becomes nonzero. 
Moreover, the solutions produced by the semidefinite 
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relaxation are physically meaningless in those cases. 
Indeed, examples of nonconvexity have long been dis- 
cussed in the literature, e.g., [21]-[23]. Hence it is 
important to develop systematic methods for solving 
OPF involving mesh networks when convex relaxation 
fails. See, e.g., [24] for branch-and-bound algorithms for 
solving OPF when the duality gap is nonzero. 

The papers above are all based on the bus injection 
model. In this paper, we introduce a branch flow model 
on which OPF and its relaxations can also be defined. 
Our model is motivated by a model first proposed by 
Baran and Wu in [25], [26] for the optimal placement and 
sizing of switched capacitors in distribution circuits for 
Volt/VAR control. As explained in Section III below, one 
of the key insights we highlight here is that the Baran- 
Wu model of [25], [26] can be treated as a particular 
relaxation of our branch flow model where the phase 
angles of the voltages and currents are ignored. By 
recasting their model as a set of linear and quadratic 
equality constraints, [27], [28] observe that relaxing the 
quadratic equality constraints to inequality constraints 
yields a second-order cone program (SOCP). It proves 
that the SOCP relaxation is exact for radial networks, 
when there are no upper bounds on the loads. This 
result is extended here to mesh networks with line limits, 
and convex, as opposed to linear, objective functions 
(Theorem 1). See also [29], [30] for various convex 
relaxations of approximations of the Baran- Wu model 
for radial networks. 

Other branch flow models have also been studied, e.g., 
in [31]-[33], all for radial networks. Indeed, [31] studies 
a similar model to that in [25], [26], using receiving- 
end branch powers as variables instead of sending-end 
branch powers as in [25], [26]. Both [32] and [33] 
eliminate voltage angles by defining real and imaginary 
parts of ViVj^ as new variables and defining bus power 
injections in terms of these new variables. This results 
in a system of linear quadratic equations in power 
injections and the new variables. While [32] develops 
a Newton-Raphson algorithm to solve the bus power 
injections, [33] solves for the branch flows through an 
SOCP relaxation for radial networks, though no proof of 
optimality is provided. 

This set of papers [25]-[33] all exploit the fact that 
power flows can be specified by a simple set of linear and 
quadratic equalities if voltage angles can be eliminated. 
Phase angles can be relaxed only for radial networks 
and generally not for mesh networks, as [34] points out 
for their branch flow model, because cycles in a mesh 
network impose nonconvex constraints on the optimiza- 
tion variables (similar to the angle recovery condition in 
our model; see Theorem 2 below). For mesh networks. 



[34] proposes a sequence of SOCP where the nononvex 
constraints are replaced by their linear approximations 
and demonstrates the effectiveness of this approach using 
seven network examples. In this paper we extend the 
Baran- Wu model from radial to mesh networks and use 
it to develop a solution strategy for OPF, as we now 
summarize. 

B. Summary 

Our purpose is to develop a formal theory of branch 
flow model for the analysis and optimization of mesh as 
well as radial networks. As an illustration, we formulate 
OPF within this alternative model, propose relaxations, 
characterize when a relaxed solution is exact, prove that 
our relaxations are always exact for radial networks when 
there are no upper bounds on loads but may not be exact 
for mesh networks, and show how to use phase shifters 
to convexify a mesh network so that a relaxed solution 
is always optimal for the convexified network. A similar 
set of results have been proved in the sequence of papers 
[14], [16]-[19] for the bus injection model, even though 
the results have distinct characters in each model and the 
proof techniques are completely different. 

Specifically, we first formulate in Section II the OPF 
problem using branch flow equations involving complex 
bus voltages and complex branch current and power 
flows. In Section III we describe our solution approach 
that consists of two relaxation steps, as illustrated in 
Figure 1: 

• Angle relaxation: relax OPF by eliminating voltage 
and current angles from the branch flow equations. 
This yields the (extended) Baran- Wu model and a 
relaxed problem OPF-ar which is still nonconvex. 

• Conic relaxation: relax OPF-ar by taking the convex 
hull of the feasible set. This yields a second-order 
cone program OPF-cr that can be solved efficiently. 

In Section IV we prove that the conic relaxation OPF- 
cr is always exact even for mesh networks, provided there 
are no upper bounds on real and reactive loads, i.e., any 
optimal solution of OPF-cr is also optimal for OPF-ar. 
Given an optimal solution of OPF-ar, whether we can 
derive an optimal solution of the original OPF depends 
on whether we can recover the voltage and current angles 
correctly from the given OPF-ar solution. In Section V 
we characterize the exact condition (the angle recovery 
condition) under which this is possible, and present two 
angle recovery algorithms. It turns out that the angle 
recovery condition always holds for a radial network 
and hence solving OPF-cr always produces an optimal 
solution for OPF. For a mesh network, the angle recovery 
condition corresponds to the requirement that the implied 
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Fig. 1: Proposed solution strategy for solving OPF. 



phase angle differences sum to zero around every loop. 
The given OPF-ar solution may not satisfy this condition, 
but our characterization can be used to check if it yields 
an optimal solution for OPF. 

These results suggest an algorithm for solving OPF as 
summarized in Figure 2. 
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radial 
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OPF solution 

Fig. 2: Proposed algorithm for solving OPF (11)-(12) 
without phase shifters. The details are explained in 
Sections II-V. 

If a relaxed solution for a mesh network does not 
satisfy the angle recovery condition, then it is infeasible 
for OPF. In Part II of this paper, we propose a simple 
way to convexify a mesh network using phase shifters 
so that any relaxed solution of OPF-ar can be mapped to 
an optimal solution of OPF for the convexified network, 
with an optimal cost that is lower than or equal to that 
of the original network. 



C. Extensions: radial networks and equivalence 

In [35], [36], we prove a variety of sufficient condi- 
tions under which the conic relaxation proposed here is 
exact for radial networks. The main difference from The- 
orem 1 below is that, [35], [36] allow upper bounds on 
the loads but relax upper bounds on voltage magnitudes. 
Unlike the proof for Theorem 1 here, those in [35], [36] 
exploit the duality theory. 

The bus injection model and the branch flow model 
are defined by different sets of equations in terms of 
their own variables. Each model is self-contained: one 
can formulate and analyze power flow problems within 
each model, using only nodal variables or only branch 
variables. Both models (i.e., the sets of equations in 
their respective variables), however, are descriptions of 
the Kirchhoff's laws. In [37] we prove formally the 
equivalence of these models, in the sense that given a 
power flow solution in one model, one can derive a 
corresponding power flow solution in the other model. 
Although the semidefinite relaxation in the bus injection 
model is very different from the convex relaxation pro- 
posed here, [37] also establishes the precise relationship 
between the relaxations of these two models. 

As a consequence, the various conditions on radial net- 
works proved in [16]-[19] for the exactness of semidef- 
inite relaxation imply immediately the exactness of the 
relaxations proposed here. Conversely, the conditions in 
[35], [36] for the conic relaxation immediately imply that 
there is a matrix solution of the semidefinite relaxation 
in [16]-[19] that is of rank 1. This is useful because 
some results are easier to formulate and prove in one 
model than in the other. For instance, it is hard to see 
how the upper bounds on voltage magnitudes and the 
technical conditions on the line impedances in [35], [36] 
for exactness in the branch flow model affect the rank 
of the semidefinite matrix variable in the bus injection 
model, although [37] guarantees their equivalence. 

II. Branch flow model 

Let R denote the set of real numbers and C denote the 
set of complex numbers. A variable without a subscript 
usually denotes a vector with appropriate components, 
e.g., s := {si,i = 1, . . . , n), 5 := {Sij, (ij) G E). For 
a complex scalar or vector a, a* denotes its complex 
conjugate. For a vector a = (ai, . . . ,a/e), a-i denotes 
(ai, . . . , a^-i, a^+i, ttjt). For a matrix A, denotes its 
transpose and A* its complex conjugate transpose. Given 
a directed graph G = (A^, E), denote a link in Eby (i^j) 
or i ^ j if it points from node i to node j. We will use 
e, (i, j), or i ^ j interchangeably to refer to a link in 
E. We write i ^ j ifi and j are connected, i.e., if either 
(i, j) ^ E ox (j, i) ^ E (but not both). 



4 



A. Branch flow model 

Let G = [N^ E) be a connected graph representing a 
power network, where each node in represents a bus 
and each Unk in E represents a Une (condition Al). We 
index the nodes by i = 0, 1, . . . , n. The power network 
is called radial if its graph G is a tree. For a distribution 
network, which is typically radial, the root of the tree 
(node 0) represents the substation bus. For a (generally 
meshed) transmission network, node represents the 
slack bus. 

We regard G as a directed graph and adopt the follow- 
ing orientation for convenience (only). Pick any spanning 
tree T := {N,Et) of G rooted at node 0, i.e., T is 
connected and Et Q E has n links. All links in Et point 
away from the root. For any link in \ Et that is not in 
the spanning tree T, pick an arbitrary direction. Denote 
a link by (i, j) or i ^ j if it points from node i to node 
j. Henceforth we will assume without loss of generality 
that G and T are directed graphs as described above. ^ For 
each link (i, j) G E, let zij = rij + \Xij be the complex 
impedance on the line, and yij := =: Qij — ibij be 
the corresponding admittance. For each node i E let 

= n + ixi be the shunt impedance from i to ground, 
and yi := 1/zi =: gi - ibi. 

For each (i, j) G E, let lij be the complex current 
from buses i to j and Sij = Pij+iQij be the sending-end 
complex power from buses i to j. For each node i G A^, 
let Vi be the complex voltage on bus i. Let Si be the 
net complex power injections, which is generation minus 
load on bus i. We use Si to denote both the complex 
number pi + iqi and the pair (pi^Qi) depending on the 
context. 

As customary, we assume that the complex voltage Vq 
is given and the corresponding complex net generation sq 
is a variable. For power flow analysis, we assume other 
power injections s := {si^i = 1, . . . ,n) are given. For 
optimal power flow, VAR control, or demand response, 
s are control variables as well. 

Given z := {zij,{i,j) G E, zi^i G N) and 
bus power injections 5, the variables (5, /, y, 5o) := 
{Sij, lij, (ij) G E, Vi,i = sq) satisfy the 

Ohm's law: 

Vi-Vj = zijiij, y{ij)eE (1) 

the definition of branch power flow: 

Sij = Vil*j, y{i,j)eE (2) 

^The orientation of G and T are different for different spanning 
trees T, but we often ignore this subtlety in this paper. 



and power balance at each bus: for all j G N, 

J2 S,k -Y.^Si, - z,,|I,,f ) + = sj (3) 

We will refer to (l)-(3) as the branch flow 
model/equations. Recall that the cardinality |A^| = n + 1 
and let \E\ =: m. The branch flow equations (l)-(3) 
specify 2m + n + 1 nonlinear equations in 2m + n + 1 
complex variables {S,I,V,sq), when other bus power 
injections s are specified. 

We will call a solution of (l)-(3) a branch flow 
solution with respect to a given 5, and denote it by 
x{s) := (5,/,y,5o). Let X{s) C C^^^+^+i be the set 
of all branch flow solutions with respect to a given s: 

X{s) := {x := {S, /, V, sq) \ x solves (l)-(3) given s} 

(4) 

and let X be the set of all branch flow solutions: 

X := U X{s) (5) 

For simplicity of exposition, we will often abuse notation 
and use X to denote either the set defined in (4) or 
that in (5), depending on the context. For instance, X 
is used to denote the set in (4) for a fixed s in Section V 
for power flow analysis, and to denote the set in (5) in 
Section IV for optimal power flow where s itself is also 
an optimization variable. Similarly for other variables 
such as X for x{s), 

B. Optimal power flow 

Consider the optimal power flow problem where, 
in addition to {S, I,V, sq), s is also an optimization 
variable. Let pi pf — Pi and qi '-^ qf — qi where p^ 
and qf (p^ and q^) are the real and reactive power gen- 
eration (consumption) at node i. For instance, [25], [26] 
formulate a VoltA^AR control problem for a distribution 
circuit where qf represent the placement and sizing of 
shunt capacitors. In addition to (l)-(3), we impose the 
following constraints on power generation: for i G N, 

f,<pl<pf, ql<q!<ql (6) 

In particular, any of pf^qi can be a fixed constant by 
specifying that p^ = pf and/or = • For instance, in 
the inverter-based VAR control problem of [27], [28], pf 
are the fixed (solar) power outputs and the reactive power 
qf are the control variables. For power consumption, we 
require, for i G iV, 

p^<pl<pt, ti<qt<q'l (7) 
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The voltage magnitudes must be maintained in a tight 
range: for i = 1, . . . , n, 

Vi < \Vi\^ < Vi (8) 

Finally, we impose line flow limits: for all (i, j) G E, 

I Sij I < Sij (9) 

We allow any objective function that is convex and 
does not depend on the angles ZVi^ /Llij of voltages and 
currents. For instance, suppose we aim to minimize real 
power losses r^j|/^jp, minimize real power generation 
costs Cip^, and maximize energy savings through con- 
servation voltage reduction (CVR). Then the objective 
function takes the form (see [27], [28]) 

for some given constants > 0. We also allow 

the cost to be quadratic in real power as is commonly 
assumed. 

To simplify notation, let iij := and Vi := 
Let s3 := {sf, i = l,...,n) = (pf,^f, i = l,...,n) be 
the power generations, and := (5^, i = 1, . . . , n) = 
{ph^h ^ ~ l?---?^) the power consumptions. Let 5 
denote either s^ — s^ or {s^ ^s^) depending on the context. 
Given a branch flow solution x := ^(5) := (5, /, sq) 
with respect to a given 5, let y := y{s) := (5, ^, s^) de- 
note the projection of x that have phase angles Zl^, Z/^^ 
eliminated. This defines a projection function h such that 
y = to which we will return in Section IIL Then 

our objective function is / [h{x)^ . We assume / (y, s) 
is convex (condition A2); in addition, we assume / is 
strictly increasing in (i, j) G E and nonincreasing in 
load (condition A3). Finally, let 

S := { {S, V, 50, s) I {v, 50, s) satisfies (6) - (9) } 

All quantities are optimization variables, except Vq 
which is given. 

The optimal power flow problem is 
OPF: 

mm f(h[x),s) (11) 

subject to G X, (S, sq, 5) G S (12) 

where X is defined in (5). To avoid triviality, we assume 
the problem is feasible (condition A4). 

The feasible set is specified by the nonlinear branch 
flow equations and hence OPF (11)-(12) is in general 
nonconvex and hard to solve. The goal of this paper is 
to propose an efficient way to solve OPF by exploiting 
the structure of the branch flow model. 



C. Notations and assumptions 

The main variables and assumptions are summarized 
in Table I and below for ease of reference: 



TABLE L Notations. 



T 


(directed) network graph G and a span- 
ning tree T of G 


B, Bt 


reduced (and transposed) incidence 
matrix of G and the submatrix corre- 
sponding to T 


Vi, Vi 


complex voltage on bus i with Vi := 


Si =Pi-\- IQi 

Pi=pl - Pi 

Qi = q! - Qi 


net complex load power on bus i 
net real power equals generation minus 
load; 

net reactive power equals generation 
minus load 




complex current from buses i to j with 

iij :— \Iij\ 


Sij — Pij -\- iQij 


complex power from buses i to j 
(sending-end) 


X 

Y 
Y 


set 01 all branch now solutions that 
satisfy (l)-(3) either for some s, or 
for a given s (sometimes denoted more 
accurately by X(s)); 
set of all relaxed branch flow solutions 
that satisfy (13)-(16) either for a given 
s or for some s; 

convex hull of Y, i.e., solutions of 
(13)-(15) and (22); 


x = (S',/,y,so) gx 
y = {SJ,v,so) G Y 
y = h{x)] x = he{y) 


vector X of power flow variables 

and its projection y; 

projection mapping y and an inverse he 


Zij , Zi 


impedance on line connecting buses i 
and j and shunt implemance from bus 
i to ground, respectively 


f 


objective function of OPF, of the form 
f(h{x),s) 



Al The network graph G is connected. 
A2 The cost function /(y, s) for optimal power flow is 
convex. 

A3 The cost function f{y^s^) is strictly increasing in 

(i, j) G E and nonincreasing in load s^. 
A4 The optimal power flow problem OPF (11)-(12) is 
feasible. 

These assumptions are standard and realistic. For in- 
stance, the objective function in (10) satisfies conditions 
A2-A3. A3 holds if the cost function is strictly increas- 
ing in power loss in the lines. 

III. Relaxations and solution strategy 
We now describe our solution approach. 

A. Relaxed branch flow model 

Substituting (2) into (1) yields Vj = Vi- ZijSlj/V-'. 
Taking the magnitude squared, we have vj = Vi + 
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\zij\'^iij - {zijSl- + zl-Sij). Using (3) and (2) and in 
terms of real variables, we therefore have 

Pj = Pjk ~ ^^^^^^^ + (^^^ 

^3 = X] -^^Qij - XiAj) + ^3^3^ Vj(14) 

Vj = livijPij -\- XijQij^ -\- {j'ij ~h 

y{i,j)eE (15) 
= — v(i,i)G£; (16) 

We will refer to (13)-(16) as the relaxed (branch 
flow) model/equations and a solution a relaxed (branch 
flow) solution. These equations were first proposed 
in [25], [26] to model radial distribution circuits. 
They define a system of equations in the variables 
{P,QJ,v,po,qo) := {Pij^QijJij, (ij) ^ E, Vi,i = 
po^qo). We often use {S^I^v^sq) as a short- 
hand for (P^Q^i^v^po^qo). The relaxed model has a 
solution under A4. 

In contrast to the original branch flow equations (1)- 
(3), the relaxed equations (13)-(16) specifies 2(m+n+l) 
equations in 3m + n + 2 real variables (P, Q, ^, ^,Po, ^o), 
for a given 5. For a radial network, i.e., G is a tree, 
m = \E\ = \N\ — 1 = n. Hence the relaxed system (13)- 
(16) specifies 4n + 2 equations in 4n + 2 real variables. 
It is shown in [38] that there are generally multiple 
solutions, but for practical networks where |Vb| — 1 and 
Hj, Xij are small p.u., the solution of (13)-(16) is unique. 
Exploiting structural properties of the Jacobian matrix, 
efficient algorithms have also been proposed in [39] to 
solve the relaxed branch flow equations. 

For a connected mesh network, m = \E\ > |A^| — 1 
n, in which case there are more variables than equations 
for the relaxed model (13)-(16), and therefore the so- 
lution is generally nonunique. Moreover, some of these 
solutions may be spurious, i.e., they do not correspond to 
a solution of the original branch flow equations (l)-(3). 

Indeed, one may consider (S^I^v^sq) as a projection 
of (5, /, y, So) where each variable lij or Vi is relaxed 
from a point in the complex plane to a circle with a 
radius equal to the distance of the point from the origin. 
It is therefore not surprising that a relaxed solution of 
(13)-(16) may not correspond to any solution of (1)- 
(3). The key is whether, given a relaxed solution, we 
can recover the angles ZVi^ZIij correctly from it. It 
is then remarkable that, when G is a tree, indeed the 
solutions of (13)-(16) coincide with those of (l)-(3). 
Moreover /or a general network, (13)-(16) together with 
the angle recovery condition in Theorem 2 below are 



indeed equivalent to (l)-(3), as explained in Remark 5 
of Section V. 

To understand the relationship between the branch 
flow model and the relaxed model and formulate our 
relaxations precisely, we need some notations. Fix an 
s. Given a vector (5,/,l/,5o) G C^^+^+i, define its 

projection h : C^^+^+l ^ ^3m+n+2 f^^g^ _ 

{P,QJ,v,po,qo) where 

Pi = Re Si, qi = Im Si, Vi = \Vi\'^ (18) 

Let Y C C^rn+n+l ^^^^^^ ^j^^ ^jj y 




Fig. 3: X is the set of branch flow solutions and 
Y = h{Y) is the set of relaxed solutions. The inverse 
projection kg is defined in Section V. 

whose projections are the relaxed solutions:^ 

Y := [y := {S,I so)\h{y) solves (13)-(16)} (19) 
Define the projection Y := /i(Y) of Y onto the space 

Y := { ^ := (5, 1, V, so) \ y solves (13)-(16) } 
Clearly 

X C Y and h{X) C h{Y) =: Y 
Their relationship is illustrated in Figure 3. 

B. Two relaxations 

Consider the OPF with angles relaxed: 

min / ( h(x), s ] 

x,s \ J 

subject to X G Y, (5, v^sq^s) G S 

Clearly, this problem provides a lower bound to the 
original OPF problem since Y ^ X. Since neither h{x) 

^As mentioned earlier, the set defined in (19) is strictly speaking 
Y(s) with respect to a fixed s. To simplify exposition, we abuse 
notation and use Y to denote both Y(s) and Usec^ ^i^)^ depending 
on the context. The same applies to Y and Y etc. 
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nor the constraints in Y involves angles ZVi^ZIij, 
problem is equivalent to the following 
OPF-ar: 



this 



mm 



subject to 



fiy.s) (20) 

yeY, {S,v,so,s)eS (21) 



The feasible set of OPF-ar is specified by a system of 
linear-quadratic equations. Hence OPF-ar is in general 
still nonconvex and hard to solve directly. The key to 
our solution is the observation that the only source of 
nonconvexity is the quadratic equalities in (16). Relax 
them to inequalities: 



> 



^ij + Qij 



{iJ)eE (22) 

:m+n+l y aS 



and define the convex hull Y C 

Y := {y := {SJ,v,so)\y solves (13)-(15) and (22)} 

Consider the following relaxation of OPF-ar: 
OPF-cr: 



mm 



subject to 



fiy.s) (23) 
yeY, {S,v,so,s)eS (24) 



Clearly OPF-cr provides a lower bound to OPF-ar since 
YD Y. 

C. Solution strategy 

In the rest of this paper, we will prove the following 
results: 

1) OFP-cr is convex. Moreover, if there are no upper 
bounds on loads, then the conic relaxation is exact 
so that any optimal solution (^cr, Scr) of OPF-cr is 
also optimal for OPF-ar for mesh as well as radial 
networks (Section IV, Theorem 1). 

2) Given a solution (^^r, Sar) of OPF-ar, if the network 
is radial, then we can always recover the phase 
angles ZVi, Zlij uniquely to obtain an optimal solu- 
tion (x*, 5*) of the original OPF (11)-(12) through 
an inverse projection (Section V, Theorems 2 and 
4). 

3) For a mesh network, an inverse projection may 
not exist to map the given (yiar^Sar) to a feasible 
solution of OPF. Our characterization can be used 
to determined if {yar^ Sar) turns out to be globally 
optimal. 

These results motivate the algorithm in Figure 2. 

In Part II of this paper, we show that a mesh network 
can be convexified so that {yar^Sar) can always be 
mapped to an optimal solution of OPF for the convex- 
ified network. Moreover, convexification requires phase 



shifters only on lines outside an arbitrary spanning tree 
of the network graph. 



IV. Exact conic relaxation 

Our first key result says that OPF-cr is exact. 

Theorem 1: Suppose — = oo, i ^ N . Then 
OPF-cr is convex. Moreover, it is exact, i.e., any optimal 
solution of OPF-cr is also optimal for OPF-ar. 

Proof: The feasible set is convex since the nonlinear 
inequalities in Y can be written as the following second 
order cone constraint: 



2Qij 



Since the objective function is convex, OPF-cr is a conic 
optimization.^ To prove that the relaxation is exact, it 
suffices to show that any optimal solution of OPF-cr 
attains equality in (22). 

Assume for the sake of contradiction that (y*, «§*) := 
5 '^>kO' '^*0' '^*' ^i) Optimal for OPF-cr, but a 
link (i, j) G E has strict inequality, i.e., > 
[P*]ij^+[(5*]^j^. For some e > to be determined below, 
consider another point (y^s) = (S^I^v^Sq^Sq^s^^s^) 
defined by: 



V = 
£ij = 
Sij = 



S—iq — 



«1 



[S:^]—ij 

+ %-^/2 



where a negative index means excluding the indexed 
element from a vector. Since £ij = [£^]ij — s) has 
a strictly smaller objective value than (y*,^*) because 
of assumption A3. If s) is a feasible point, then it 
contradicts the optimality of 

It suffices then to check that there exists an e > such 
that {y, s) satisfies (6)-(9), (13)-(15) and (22), and hence 
is indeed a feasible point. Since (^*, s^) is feasible, (6)- 
(8) hold for s) too. Similarly, s) satisfies (13)-(14) 
at all nodes k ^ i^j and (9), (15), (22) across all links 
(A;,/) ^ [hj)- We now show that satisfies (13)- 

(14) also at nodes i, j, and (9), (15), (22) across link 

{hi)- 

^If the objective function is linear, such as (10), then OPF-cr is 
an SOCP. This is the case proved in [27] for radial networks. This 
result is extended here to mesh networks with line limits and convex 
objective functions. 
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Proving (13)-(14) is equivalent to proving (3). At 
node i, we have 

h = sf-st = [si]i-[s';]i-Zije/2 

+yiVi - Zij£/2 
— ^ ^ij' + (^Sij + Zije/2^ 



k^i 



i^y k^i 

At node j, we have 

% = ~S]-~S) = [si]j-[sl]j-Zij£/2 

j^k i'^j 

+y^Vj - Zijs/2 
— 5^ ^jk — ^ (^Si'j — Zi'jii'j^ + yjVj 

j^k i'^j^i'^i 

{Sij + Zije/2) - Zij{iij + e)^ - Zije/2 
= ^ Sjk — ^ [Si'j — Zi'jii'j^ + yjVj 

j^k i'^j 

Hence (13)-(14) hold at nodes i, j. 

For (9) across link (i, j), we have \Sij\^ = \Sij\^ — 

I (2(4.5,, + z,,S*^ - £|z,,f ) < |5,,f < 4. for 

small enough e > 0. 

For (15) across link (i, j), we have 



Vj — [l^-K ] 2 2 ( Tij ] 2 ji H~ Xij \Q * ] 2 ji ) 



• For (22) across link (i, j), we have 

v-l- -P^ -O^ 
= [v*]i {[Qij -£)- i[P*]ij - rije/2f 

= {[vM^^J - [PAi - ml) 

+ £(4+4)/4) 

Since [v^]i[i^]ij-[P^]fj-[Q^]fj > 0, we can choose 
an £ > sufficiently small such that iij > [Pf^ + 

This completes the proof. ■ 



Remark 1: Assumption A3 is used in the proof here 
to contradict the optimality of (^*, 5*). Instead of A3, if 
f{y^s) is nondecreasing in i, the same argument shows 
that, given an optimal with a strict inequality 

> [P^]ij'^ + one can choose ^ > to 

obtain another optimal point s) that attains equality 
and has a cost /(^, 5) < /(^*, s*). Hence, in the absence 
of A3, there is always an optimal solution of OPF-cr that 
is also optimal for OPF-ar, even though it is possible 
that the convex relaxation OPF-cr may also have other 
optimal points with strict inequality that are infeasible 
for OPF-ar. 

Remark 2: The condition in Theorem 1 is equivalent 
to the "over-satisfaction of load" condition in [14], [17]. 
It is needed because we have increased the loads on 
buses i and j to obtain the alternative feasible solution 
(y, 5). As we show in the simulations in [40], it is 
sufficient but not necessary. See also [35], [36] for 
exact conic relaxation of OPF-cr for radial networks 
where this condition is replaced by an alternative set 
of assumptions. 

V. Angle relaxation 

Theorem 1 justifies solving the convex problem OPF- 
cr for an optimal solution of OPF-ar. Given a solution 
{y, s) of OPF-ar, when and how can we recover a 
solution (x, 5) of the original OPF (11)-(12)? It depends 
on whether we can recover a solution x to the branch 
flow equations (l)-(3) from y, given any s. 

Hence, for the rest of Section V, we fix an 5. We abuse 
notation in this section and write x, ^, 6^ X, Y, Y instead 
of , ^(5) , e{s) , X{s) , Y(5) , Y(5) respectively. 

A. Angle recovery condition 

Fix a relaxed solution y := (S^I^v^sq) G Y. Define 
the (n + 1) X m incidence matrix C of G by 



1 if link e leaves node i 
Cie = < — 1 if link e enters node i 
otherwise 



(25) 



The first row of C corresponds to node 0, the refer- 
ence bus with a given Vq = |Vo|e*^°- In this paper 
we will only work with the m x n reduced incidence 
matrix B obtained from C by removing the first row 
(corresponding to Vq) and taking the transpose, i.e., for 
e G E', i = 1, . . . , n. 



Bei = <-l 



if link e leaves node i 
if link e enters node i , 
otherwise 
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Since G is connected, m > n and rank(5) = n [41]. 
Fix any spanning tree T = (A^, Et) of G. We can 
assume without loss of generality (possibly after re- 
labeling some of the links) that Et consists of links 
e = 1, . . . , n. Then B can be partitioned into 



B 



Bj" 
B± 



(26) 



where the n x n submatrix Bt corresponds to hnks in T 
and the {m — n)xn submatrix B± corresponds to Unks 

in :=G\ T. 

Let (3 := I3{y) G [-7r,7r]"^ be defined by: 



^3 



Z{vi-z*jSij), ii,j)eE (27) 



Informally, Pij is the phase angle difference across link 
that is impUed by the relaxed solution y. Write /? 

as 



(28) 



where /3t is n x 1 and I3± is {m — n) x 1. 

Recall the projection mapping h : C^^+^+^ ^ 
]^3m+n+2 defined in (17)-(18). For each 9 := (6>i,i = 
l,...,n) G [— 7r,7r]^, define the inverse projection 

he : R3m+n+2 _^ (C2m+n+l h0{P, Q J, V , po, Qo) = 

(5, /, y, 5o) where 



Sij 




(29) 


Uj 




(30) 


Vi 


:= v^e''^ 


(31) 




:= po + i^o 


(32) 



These mappings are illustrated in Figure 3. 

By definition of /i(X) and Y, a branch flow solution 
in X can be recovered if a given relaxed solution y is in 
/i(X) and cannot be recovered if y is in Y\/i(X). In other 
words, /z(X) consists of exactly those points y G Y for 
which there exist Q such that their inverse projections 
heiji) are in X. Our next key result characterizes the 
exact condition under which such an inverse projection 
exists, and provides an explicit expression for recovering 
the phase angles ZT^, Z/^^ from the given y. 

A cycle c in G is an ordered list c = {ii^ . . . ^ik) 
of nodes in N such that (ii ^ 22), . . . , (ik ^ h) are all 
links in E. We will use '{i^j) G c' to denote a link i ~ j 
in the cycle c. Each such link i ^ j may be in the same 
orientation ((i, j) G or in the opposite orientation 
((j,i) G Let ^ be the extension of /3 from directed 
links to undirected links: if (i, j) G £^ then := f3ij 
and := — Aj. For any d-dimensional vector a, let 
V{a) denote its projection onto [— 7r,7r]^. 



Theorem 2: Let T be <2ny spanning tree of G. Con- 
sider a relaxed solution ^ G Y and the corresponding ^ 
defined by (27)-(28) in terms of y. 
1) There exists a unique 0* G [— 7r,7r]^ such that 
(y) is a branch flow solution in X if and only if 



B^B-^Pt = P± 



(33) 



2) The angle recovery condition (33) holds if and only 
if for every cycle c in G 



E 



= 



(34) 



3) If (33) holds then O^^V {B-^^t)- 

Remark 3: Given a relaxed solution y := (5, ^, i;, 5o), 
Theorem 2 prescribes a way to check if a branch 
flow solution can be recovered from it, and if so, the 
required computation. The angle recovery condition (33) 
is a condition on y and depends only on the network 
topology through the reduced incidence matrix B. The 
choice of spanning tree T corresponds to choosing n 
linearly independent rows of B to form Bt and does 
not affect the conclusion of the theorem. 

Remark 4: When it holds, the angle recovery condi- 
tion (34) has a familiar interpretation (due to Lemma 3 
below): the voltage angle differences (implied by y) sum 
to zero around any cycle. 

Remark 5: A direct consequence of Theorem 2 is that 
the relaxed branch flow model (13)-(16) together with 
the angle recovery condition (33) is equivalent to the 
original branch flow model (l)-(3). That is, x satisfies 
(l)-(3) if and only if y = h{x) satisfies (13)-(16) and 
(33). The challenge in computing a branch flow solution 
X is that (33) is nonconvex. 

The proof of Theorem 2 relies on the following 
important lemma that gives a necessary and sufficient 
condition for an inverse projection ho{y) defined by 
(29)-(32) to be a branch flow solution in X. 

Lemma 3: Given any y := (5, sq) in Y: 

1) There exists at most one 9 that satisfies 

Oi-Oj = Z{vi-z:jSij), {iJ)eE (35) 

2) Let 9 be the unique solution of (35). Then h^{y) is 
a branch flow solution in X if 

9i = 9i + 2ki7T for all A;^ G {±0,±1,...,} 

Conversely if h^{y) is in X then 

9i = 9i + 2ki7r for some ki G {±0, ±1, . . . , } 

Proof: In vector form, the condition (35) is equiv- 
alent to: whether there exist 9 that solves 



B9 = p 



(36) 
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i.e., whether /? is in the range space of B. Since G is 
connected, m> n and rank(A) = n and hence (36) has 
at most one solution. 

Suppose 6 is the unique solution of (35). We need 
to show that (13)-(16) together with (29)-(32) and (35) 
implies (l)-(3). Now (13) and (14) are equivalent to (3). 
Moreover (16) and (29)-(31) imply (2). To prove (1), 
substitute (2) into (35) to get 

-Oj = z {v, - z^jvatj) = ^ ^^ {v^ - ^ijhjT 

Hence 

Z.Vj = Oj = Z{Vi-ZijIij) (37) 
From (15) and (2), we have 

{Vjl'^ = + \zij\'^\Iij\'^ — {zijS^j + z^jSij) 

~ I "^i Zij lij I 

This together with (37) implies Vj^Vi — Zijiij which 
is (1). Finally, it is easily verified that 9 given in the 
theorem also defines an inverse projection that maps 
y to a branch flow solution h^{y) in X. 

Suppose h^{y) E X. From (1) and (2), we have 

ViV; = -4.5,,. Then 9,-9j ^z[v,- 45,,) ± 

2/c^j7r for some integer ki. Hence 9i = 9i±2ki7r for some 
integer ki where 9 is the unique solution of (35). ■ 

Proof of Theorem 2: Since m > n and rank(fi) = 
n, we can always find n linearly independent rows of B 
to form a basis. The choice of this basis corresponds to 
choosing a spanning tree of G, which always exists since 
G is connected [42, Chapter 5]. Assume without loss of 
generality that the first n rows is such a basis so that B 
and P are partitioned as in (26) and (28) respectively. 
Then Lemma 3 implies that there exists a unique G 
[— TT, tt]^ such that ho^ (y) is a branch flow solution in X 
if and only if 9^ = V{9) where 9 solves (36), i.e.. 



Bj" 


9 = 








A. 



(38) 



Since T is a spanning tree, the n x n submatrix Bt 
has a full rank and hence is invertible. Therefore 9 = 
B^^I^T' Moreover (38) has a unique solution if and only 
if B^9 = B^B-^Pt = 

We are left to prove the equivalence of (33) and (34). 
Recall that the spanning tree T defines the orientation 
of all links in T to be directed away from the root node 
0. Let T{i ^ j) denote the unique path from node i to 
node j in T; in particular, r(0 ^ j) consists of links 
all with the same orientation as the path and T{j ^ 0) 



of links all with the opposite orientation. Then it can be 
verified directly that the inverse Bji is given by 

. _ f-1 if link e is in r(0 ^ i) 

^ ^ '~ |o otherwise 

Then B^^^t represents the (negative of the) sum of 
angle differences on the path r(0 ^ i) for each node 

iGT: 

Hence B±B^^/3t is the sum of voltage angle differences 
from node i to node j along the unique path in T, for 
every link (i, j) G E\ Et not in the tree T. To see this, 
we have, for each link e := (i, j) E E \ Et, 

[B±B-'Pt], = [B-'/3t],-[B-'/3t]. 
Since 

E i^Tl' = - E [He' 

the angle recovery condition (33) is equivalent to 

E [0Tl' + [^±h+ E [He' 

e'GT(O^z) e'GT(j^O) 

= E ^e' = 

where c(i, j) denotes the unique basis cycle (with respect 
to T) associated with each link (i,j) not in T [42, 
Chapter 5]. Hence (33) is equivalent to (34) on all basis 
cycles, and therefore it is equivalent to (34) on all cycles. 



B. Angle recovery algorithms 

Theorem 2 suggests a centralized method to compute 
a branch flow solution from a relaxed solution. 
Algorithm 1: centralized angle recovery. Given a 
relaxed solution y G Y, 

1) Choose any n basis rows of B and form Bt and 

2) Compute /3 from y and check if B^Bj}Pt = 

3) If not, then y ^ h{X)\ stop. 

4) Otherwise, compute 9^ = B^^^t- 

5) Compute hoSv) ^ ^ through (29)-(32). 
Theorem 2 guarantees that Hq^ {y), if exists, is the unique 
branch flow solution of (l)-(3) whose projection is y. 

The relations (2) and (35) motivate an alternative 
iterative procedure to compute the angles Z/^j, ZVi, and 
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hence a branch flow solution. This procedure is more 
amenable to a distributed implementation. 
Algorithm 2: distributed angle recovery. Given a 
relaxed solution y G Y, 

1) Choose any spanning tree T of G rooted at node 0. 

2) For j = 0, 1, . . . , n (i.e., as j ranges over the tree T, 
starting from the root and in the order of breadth- 
first search), for all children k with j — )► set 

Zljk := ^Vj-ZSjk (39) 
ZVk := ZVj - Z{vj - z^^Sjk) (40) 

3) For each link (j, /c) G E\ Et not in the spanning 
tree, node j is an additional parent of k in addition 
to fc's parent in the spanning tree from which ZV^ 
has already been computed in Step 2. 

a) Compute current angle using (39). 

b) Compute a new voltage angle O^, using the new 
parent j and (40). If ^ then angle 
recovery has failed; stop. 

If the angle recovery procedure succeeds in Step 3, then 
y together with these angles ZVk^^Ijk are indeed a 
branch flow solution. Otherwise, a link (j, k) not in 
the tree T has been identified where condition (34) is 
violated over the unique basis cycle (with respect to T) 
associated with link (j, k). 

C. Radial networks 

Recall that all relaxed solutions in Y \ /z(X) are 
spurious. Our next key result shows that, for radial 
network, /i(X) = Y and hence angle relaxation is always 
exact in the sense that there is always a unique inverse 
projection that maps any relaxed solution ^ to a branch 
flow solution in X (even though X / Y). 

Theorem 4: Suppose G = T is a tree. Then 

1) h{X) = Y. 

2) given any y, 9^ := V [B~^(3) always exists and is 
the unique vector in [— tt, tt]^ such that ho^ {y) G X. 
Proof: When G = T is a tree, m = n and hence 

B = Bt and /? = /3t- Moreover B is n x n and of 
fuU rank. Therefore (9* = 'P ^ [-7r,7r]^ always 

exists and, by Theorem 2, ho^{y) is the unique branch 
flow solution in X whose projection is y. Since this holds 
for any arbitrary y G Y, Y = /i(X). ■ 

A direct consequence of Theorem 1 and Theorem 4 
is that, for a radial network, OPF is equivalent to the 
convex problem OPF-cr in the sense that we can obtain 
an optimal solution of one problem from that of the other. 

Corollary 5: Suppose G is a tree. Given any optimal 
solution (y*,'S*) of OPF-cr, there exists a unique 6^ G 



[— 7r,7r]^ such that {hQ^{y^)^s^) is an optimal solution 
of OPF. 

VI. Conclusion 

We have presented a branch flow model and demon- 
strated how it can be used for the analysis and opti- 
mization of mesh as well as radial networks. Our results 
confirm that radial networks are computationally much 
simpler than mesh networks and we should exploit this 
advantage whenever we can. 

Specifically, we have proposed a solution strategy for 
OPF that consists of two steps: 

1) Compute a relaxed solution of OPF-ar by solving 
its conic relaxation OPF-cr. 

2) Recover from a relaxed solution an optimal solu- 
tion of the original OPF using an angle recovery 
algorithm, if possible. 

We have proved that, for radial networks, both steps 
are always exact, provided there are no upper bounds 
on loads, so this strategy guarantees a globally optimal 
solution. For mesh networks the angle recovery condition 
may not hold but can be used to check if a given relaxed 
solution is globally optimal. 

The branch flow model is an alternative to the bus 
injection model. It adds to our understanding of power 
flow solutions and their relaxations, especially with the 
clarification of the equivalence of these two models [37]. 
For instance. Theorem 2 implies that the number of 
power flow solutions depends only on the magnitude 
of voltages and currents, not on their phase angles. 
The branch flow model also has the advantage that its 
variables correspond directly to physical quantities, such 
as branch power and current flows, and therefore are 
often more intuitive than a semidefinite matrix in the 
bus injection model. 

Acknowledgment 

We are grateful to Mani Chandy and Lingwen Gan 
of Caltech, and Christopher Clarke, Michael Montoya, 
and Robert Sherick of the Southern California Edison 
for helpful discussions. We acknowledge the support of 
the Resnick Institute at Caltech, National Science Foun- 
dation through NetSE grant CNS 0911041, Department 
of Energy's ARPA-E through grant de-ar0000226, the 
National Science Council of Taiwan (R. O. C.) through 
grant NSC 101-31 13-P-008-001, Southern California 
Edison, Cisco, and the Okawa Foundation. 

References 

[1] Masoud Farivar and Steven H. Low. Branch flow model: 
relaxations and convexification. In 51st IEEE Conference on 
Decision and Control, December 2012. 



12 



[2] J. Carpentier. Contribution to the economic dispatch problem. 
Bulletin de la Societe Francoise des Electriciens, 3(8):43 1-447, 
1962. In French. 

[3] J. A. Momoh. Electric Power System Applications of Optimiza- 
tion. Power Engineering. Markel Dekker Inc.: New York, USA, 
2001. 

[4] M. Huneault and F. D. Galiana. A survey of the optimal power 
flow literature. IEEE Trans, on Power Systems, 6(2):762-770, 
1991. 

[5] J. A. Momoh, M. E. El-Hawary, and R. Adapa. A review of 

selected optimal power flow literature to 1993. Part I: Nonlinear 

and quadratic programming approaches. IEEE Trans, on Power 

Systems, 14(1):96-104, 1999. 
[6] J. A. Momoh, M. E. El-Hawary, and R. Adapa. A review of 

selected optimal power flow literature to 1993. Part II: Newton, 

linear programming and interior point methods. IEEE Trans. 

on Power Systems, 14(1): 105 - 111, 1999. 
[7] K. S. Pandya and S. K. Joshi. A survey of optimal power flow 

methods. /. of Theoretical and Applied Information Technology, 

4(5):450-458, 2008. 
[8] B Stott and O. Alsag. Fast decoupled load flow. IEEE Trans. 

on Power Apparatus and Systems, PAS-93(3):859-869, 1974. 
[9] O. Alsa9, J Bright, M Prais, and B Stott. Further developments 

in LP-based optimal power flow. IEEE Trans, on Power 

Systems, 5(3):697-711, 1990. 
[10] K. Purchala, L. Meeus, D. Van Dommelen, and R. Belmans. 

Usefulness of DC power flow for active power flow analysis. In 

Proc. of IEEE PES General Meeting, pages 2457-2462. IEEE, 

2005. 

[11] B. Stott, J. Jardim, and O. Alsa?. DC Power Flow Revisited. 
IEEE Trans, on Power Systems, 24(3): 1290-1300, Aug 2009. 

[12] X. Bai, H. Wei, K. Fujisawa, and Y. Wang. Semidefinite 
programming for optimal power flow problems. Infl J. of 
Electrical Power & Energy Systems, 30(6-7):383-392, 2008. 

[13] X. Bai and H. Wei. Semi-definite programming-based method 
for security-constrained unit commitment with operational and 
optimal power flow constraints. Generation, Transmission & 
Distribution, lET, 3(2): 182-197, 2009. 

[14] J. Lavaei and S. Low. Zero duality gap in optimal power flow 
problem. IEEE Trans, on Power Systems, To Appear, 27(1):92- 
107, February 2012. 

[15] J. Lavaei. Zero duality gap for classical OPF problem con- 
vexifies fundamental nonlinear power problems. In Proc. of the 
American Control Conf, 2011. 

[16] S. Sojoudi and J. Lavaei. Physics of power networks makes 
hard optimization problems easy to solve. In IEEE Power & 
Energy Society (PES) General Meeting, July 2012. 

[17] S. Bose, D. Gayme, S. H. Low, and K. M. Chandy. Optimal 
power flow over tree networks. In Proc. Allerton Conf. on 
Comm., Ctrl, and Computing, October 2011. 

[18] B. Zhang and D. Tse. Geometry of feasible injection region of 
power networks. In Proc. Allerton Conf. on Comm., Ctrl, and 
Computing, October 2011. 

[19] S. Bose, D. Gayme, S. H. Low, and K. M. Chandy. Quadrat- 
ically constrained quadratic programs on acyclic graphs with 
application to power flow. arXiv:1203.5599vl, March 2012. 

[20] B. Lesieutre, D. Molzahn, A. Borden, and C. L. DeMarco. 
Examining the limits of the application of semidefinite program- 
ming to power flow problems. In Proc. Allerton Conference, 
2011. 

[21] I. A. Hiskens and R. Davy. Exploring the power flow solution 
space boundary. IEEE Trans. Power Systems, 16(3):389-395, 
2001. 

[22] B. C. Lesieutre and I. A. Hiskens. Convexity of the set of 
feasible injections and revenue adequacy in FTR markets. IEEE 
Trans. Power Systems, 20(4): 1790-1798, 2005. 



[23] Yuri V. Makarov, Zhao Yang Dong, and David J. Hill. On 

convexity of power flow feasibility boundary. IEEE Trans. 

Power Systems, 23(2):81 1-813, May 2008. 
[24] Dzung T. Phan. Lagrangian duality and branch- and-bound 

algorithms for optimal power flow. Operations Research, 

60(2):275-285, March/April 2012. 
[25] M. E. Baran and F. F Wu. Optimal Capacitor Placement 

on radial distribution systems. IEEE Trans. Power Delivery, 

4(l):725-734, 1989. 
[26] M. E Baran and F. F Wu. Optimal Sizing of Capacitors Placed 

on A Radial Distribution System. IEEE Trans. Power Delivery, 

4(l):735-743, 1989. 
[27] Masoud Farivar, Christopher R. Clarke, Steven H. Low, and 

K. Mani Chandy. Inverter var control for distribution systems 

with renewables. In Proceedings of IEEE SmartGridComm 

Conference, October 2011. 
[28] Masoud Farivar, Russell Neal, Christopher Clarke, and 

Steven H. Low. Optimal inverter VAR control in distribution 

systems with high pv penetration. In IEEE Power & Energy 

Society (PES) General Meeting, July 2012. 
[29] Joshua Adam Taylor. Conic Optimization of Electric Power 

Systems. PhD thesis, MIT, June 2011. 
[30] Joshua A. Taylor and Franz S. Hover. Convex models of dis- 
tribution system reconfiguration. IEEE Trans. Power Systems, 

2012. 

[31] R. Cespedes. New method for the analysis of distribution 
networks. IEEE Trans. Power Del, 5(l):391-396, January 
1990. 

[32] A. G. Exposito and E. R. Ramos. Reliable load flow technique 
for radial distribution networks. IEEE Trans. Power Syst., 
14(13): 1063-1069, August 1999. 

[33] R.A. Jabr. Radial Distribution Load Flow Using Conic Pro- 
gramming. IEEE Trans, on Power Systems, 21(3): 1458-1459, 
Aug 2006. 

[34] R. A. Jabr. A Conic Quadratic Format for the Load Flow 
Equations of Meshed Networks. IEEE Trans, on Power Systems, 
22(4):2285-2286, Nov 2007. 

[35] Lingwen Gan, Na Li, Ufuk Topcu, and Steven Low. Branch 
flow model for radial networks: convex relaxation. In 5 1st IEEE 
Conference on Decision and Control, December 2012. 

[36] Na Li, Lijun Chen, and Steven Low. Exact convex relaxation 
for radial networks using branch flow models. In IEEE Interna- 
tional Conference on Smart Grid Communications, November 
2012. 

[37] Subhonmesh Bose, Steven H. Low, and Mani Chandy. Equiv- 
alence of branch flow and bus injection models. In 50th 
Annual Allerton Conference on Communication, Control, and 
Computing, October 2012. 

[38] H-D. Chiang and M. E. Baran. On the existence and uniqueness 
of load flow solution for radial distribution power networks. 
IEEE Trans. Circuits and Systems, 37(3):410-416, March 1990. 

[39] Hsiao-Dong Chiang. A decoupled load flow method for distri- 
bution power networks: algorithms, analysis and convergence 
study. International Journal Electrical Power Energy Systems, 
13(3):130-138, June 1991. 

[40] Masoud Farivar and Steven H. Low. Branch flow model: 
relaxations and convexification (part II). submitted to TPS, April 
2012. 

[41] L. R. Foulds. Graph Theory Applications. Springer- Verlag, 
1992. 

[42] Norman Biggs. Algebraic graph theory. Cambridge University 
Press, 1993. Cambridge Mathematical Library. 



Branch Flow Model: Relaxations and 
Convexification (Partii) 

Masoud Farivar Steven H. Low 

Engineering and Applied Science 
Caltech 



Abstract — We propose a branch flow model for the anal- 
ysis and optimization of mesh as well as radial networks. 
The model leads to a new approach to solving optimal 
power flow (OPF) that consists of two relaxation steps. 
The first step eliminates the voltage and current angles 
and the second step approximates the resulting problem 
by a conic program that can be solved efficiently. For radial 
networks, we prove that both relaxation steps are always 
exact, provided there are no upper bounds on loads. For 
mesh networks, the conic relaxation is always exact but the 
angle relaxation may not be exact, and we provide a simple 
way to determine if a relaxed solution is globally optimal. 
We propose convexification of mesh networks using phase 
shifters so that OPF for the convexified network can always 
be solved efficiently for an optimal solution. We prove 
that convexification requires phase shifters only outside a 
spanning tree of the network and their placement depends 
only on network topology, not on power flows, generation, 
loads, or operating constraints. Part I introduces our 
branch flow model, explains the two relaxation steps, and 
proves the conditions for exact relaxation. Part II describes 
convexification of mesh networks, and presents simulation 
results. 

I. Introduction 

In Part I of this two-part paper [2], we introduce 
a branch flow model that focuses on branch variables 
instead of nodal variables. We formulate optimal power 
flow (OPF) within the branch flow model and propose 
two relaxation steps. The first step eliminates phase 
angles of voltages and currents; we call the resulting 
problem OPF-ar which is still nonconvex. The second 
step relaxes the feasible set of OPF-ar to its convex hull 
which is a (rotated) second-order cone; we call the result- 
ing problem OPF-cr. We prove that the conic relaxation 
OPF-cr is always exact even for mesh networks, provided 
there are no upper bounds on real and reactive loads, 
i.e., any optimal solution of OPF-cr is also optimal for 
OPF-ar. Given an optimal solution of OPF-ar, whether 
we can derive an optimal solution to the original OPF 
depends on whether we can recover the voltage and 

A preliminary and abridged version will appear in [1]. 



current angles correctly from the given OPF-ar solution. 
We characterize the exact condition (the angle recovery 
condition) under which this is possible, and present two 
angle recovery algorithms. It turns out that the angle 
recovery condition always holds for a radial network 
and hence solving the conic relaxation OPF-cr always 
produces an optimal solution for the original OPF. For 
a mesh network, the angle recovery condition may not 
hold, and our characterization can be used to check if a 
relaxed solution yields an optimal solution for OPF. 

In this paper, we prove that, by placing phase shifters 
on some of the branches, any relaxed solution of OPF- 
ar can be mapped to an optimal solution of OPF for 
the convexified network, with an optimal cost that is no 
higher than that of the original network. Phase shifters 
thus convert an NP-hard problem into a simpler problem. 
Our result implies that when the angle recovery condition 
holds for a relaxed branch flow solution, not only is 
the solution optimal for the OPF without phase shifters, 
but the addition of phase shifters cannot further reduce 
the cost. On the other hand, when the angle recovery 
condition is violated, then the convexified network may 
have a strictly lower optimal cost. Moreover, this benefit 
can be attained by placing phase shifters only outside an 
arbitrary spanning tree of the network graph. 

These results suggest an algorithm for solving OPF 
when there are phase shifters in mesh networks, as 
summarized in Figure 1. 

Since power networks in practice are very sparse, the 
number of lines not in a spanning tree can be relatively 
small, as demonstrated in simulations in Section V using 
the IEEE test systems with 14, 30, 57, 118 and 300 
buses, as well as a 39-bus model of a New England 
power system and two models of a Polish power system 
with more than 2,000 buses. Moreover, the placement of 
these phase shifters depends only on network topology, 
but not on power flows, generations, loads, or operating 
constraints. Therefore only one-time deployment cost 
is required to achieve subsequent simplicity in network 
operation. 

The paper is organized as follows. In Section II, we 
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Fig. 1: Proposed algorithm for solving OPF with phase 
shifters in mesh networks. The details are explained in 
this two-part paper. 

extend the branch flow model of [2] to include phase 
shifters. In Section III, we describe methods to compute 
phase shifter angles to map any relaxed solution to an 
branch flow solution. In Section IV, we explain how to 
use phase shifters to simplify OPF. In Section V, we 
present our simulation results. 



we fix an s. All quantities, such as x, X, Y, X, X^, 
are with respect to the given 5, even though that is not 
explicit in the notation. In the next section, s is also an 
optimization variable and the sets X,Y, are for 

any s. 

It is proved (Theorem 2 of [2]) that a given relaxed 
solution y can be mapped to a branch flow solution in 
X if and only if there exists an 9 that solves 

Be = ^ (6) 

where /3 := is determined by the relaxed solution: 

^3^J := Z{v,-z:^S,j), {iJ)eE (7) 

Moreover (6) has a solution and it is unique if and only 
if 

B^B-^Pt = I3± (8) 

which is equivalent to the requirement that the (implied) 
voltage angle differences sum to zero around any cycle 
c: 

where ^ij = ^ij if (ij) G E and ^ij = -^ji if (j,i) E 
E. 



II. Branch flow model with phase shifters 

We adopt the same notations and assumptions A1-A4 
of [2]. 

A. Review: model without phase shifters 

For ease of reference, we reproduce the branch flow 
model of [2] here: 

Iii^yij{Vi-Vj) (1) 
Sij = Vil*j (2) 

Recall the set X(5) of branch flow solutions given s 
defined in [2]: 

X(5) := {x := (5, /, y, 5o) I X solves (l)-(3) given s} 

(4) 

and the set X of all branch flow solutions: 

X := U X(5) (5) 

To simplify notation, we often use X to denote the set 
defined either in (4) or in (5), depending on context. In 
this section we study power flow solutions and hence 



B. Model with phase shifters 

Phase shifters can be traditional transformers or 
FACTS (Flexible AC Transmission Systems) devices. 
They can increase transmission capacity and improve 
stability and power quality [3], [4]. In this paper, we 
consider an idealized phase shifter that only shifts the 
phase angles of the sending-end voltage and current 
across a line, and has no impedance nor limits on the 
shifted angles. Specifically, consider an idealized phase 
shifter parametrized by across line (i, j), as shown 
in Figure 2. As before, let Vi denote the sending-end 




Fig. 2: Model of a phase shifter in line (i, j). 

voltage. Define Iij to be the sending-end current leaving 
node i towards node j. Let k be the point between 
the phase shifter and line impedance zij. Let Vk 
and Ik be the voltage at k and current from k to j 
respectively. Then the effect of the idealized phase shifter 
is summarized by the following modeling assumption: 

Vk = e^^- and h = hj e^^- 
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The power transferred from nodes i to j is still (defined 
to be) Sij := T^/*^ which, as expected, is equal to the 
power Vkl^ from nodes k to j since the phase shifter is 
assumed to be lossless. Applying Ohm's law across Zij, 
we define the branch flow model with phase shifters as 
the following set of equations: 



Iij = yij{Vi-Vje-"t'^^) 



Sij — Vil*j 



(9) 
(10) 

Without phase shifters = 0), (9)-(ll) reduces to the 
branch flow model (l)-(3). 

The inclusion of phase shifters modifies the network 
and enlargers the solution set of the (new) branch flow 
equations. Formally, let 

X := {x\x solves (9)-(ll) for some (/)} (12) 

All angles should be interpreted as being modulo (being 
projected into) [— 7r,7r]. Hence (/) G [— 7r,7r]^. For any 
spanning tree T of G, let "0 G T-^" stands for "0^^ = 
for all (z, j) G T", i.e., involves only phase shifters in 
branches not in the spanning tree T. Define 



Xt := { 



X I X solves (9)-(ll) for some 



GT^} 



(13) 



Since (9)-(ll) reduce to the branch flow model when 

= 0, X C Xr C X. 

III. Phase angle setting 
A. Computing cj) 

Given a relaxed solution y, there are in general many 
ways to choose angles (j) on the phase shifters to recover 
a feasible branch flow solution a; G X from y. They 
depend on the number and location of the phase shifters. 

For a network with phase shifters, we have from (9) 
and (10) 

v; 



Sij — 



Vi 



leading to F^IA* e*^^^ = - ^Ij^iy Hence Qi - 9j = 
— This changes the angle recovery condition in 
Theorem 2 of [2] from whether there exists that solves 
(6) to whether there exists {9^ (/)) that solves 



B9 = p 



(14) 



The case without phase shifters corresponds to setting 
= 0. 

We now describe two ways to compute the first 
minimizes the required number of phase shifters, and 
the second minimizes the size of phase angles. 



1 ) Minimize number of phase shifters: Our first key 
result implies that, given a relaxed solution y := 
(S^i^v^so) G Y, we can always recover a branch flow 
solution X := (5, /, sq) G X of the convexified 
network. Moreover it suffices to use phase shifters in 
branches only outside a spanning tree. This method 
requires the smallest number (m — n) of phase shifters. 

Given any (i-dimensional vector a, let V{a) denote 
its projection onto [— 7r,7r]^. 

Theorem 1: Let T be any spanning tree of G. Con- 
sider a relaxed solution ^ G Y and the corresponding /? 
defined by (7) in terms of y. 

1) There exists a unique (^*,(/)*) G [— 7r,7r]^+"^ with 
G such that ho^{y) G X^, i.e., ho^{y) is 
a branch flow solution of the convexified network. 
Specifically 

0, = V{B-^Pt) 



= V 



2) Y = X = Xt and hence Y = h{X) = hiXr). 
Prooft 

1) Write (/) = [0^ 0*^]* and set = 0. Then (14) 
becomes 



Bj" 


e = 






' ' 













The same argument as in the proof of Theorem 2 
in [2] shows that a solution (^, 0) to (14), with G 
T-^, exists and is unique if and only if 

i.e., iff 0^ = - B^B-^pT- _ 
Moreover Theorem 2 in [2] implies /^^(y) G X^. If 
(^*,0*) satisfies 



± 2fc^7r, fc^ = 0, 1, . . . 

kij = 0, 1, . . . 



± 2kij7T, 



then it is easy to verify using (9)-(ll), (13) that 
ho^{y) is also a branch flow solution in X^ with 
0* G T-^. This implies that we can always project 



into 



-TT, TT 



\n-\-m 



to obtain 



the unique vector {9 
a unique vector (^*, 0*) stated in the theorem. 
2) This follows from assertion 1. 

■ 

Remark 1: The choice of the spanning tree T does not 
affect the conclusion of the theorem. Different choices 
of T corresponds to different choices of n linearly 
independent rows of the reduced incidence matrix B 
and the resulting decomposition of B and /? into Bt 
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and Therefore T determines the phase angles 9^ 
and according to the formulae in the theorem. Since 
the objective f(h{x)^g) of OPF is independent of the 
phase angles 0*, for the same relaxed solution y, OPF- 
ps achieves the same objective value regardless of the 
choice of T. 

2) Minimize phase angles: The choice of {0^,(j)^) 
in Theorem 1 has the advantage that it requires the 
minimum number of phase shifters (only on links outside 
an arbitrary spanning tree T). It might however, require 
relatively large angles [(/)*]e at some links e outside T. 
On the other extreme, suppose we have phase shifters on 
every link. Then one can choose (6>*,(/)*) such that the 
phase shifter angles (p are minimized. 

We will formulate the problem under the simplifying 
assumption that the given relaxed solution y of OPF- 
ar and its corresponding /3 := /3{y) are such that all 
solutions {9,(f)) of (14) satisfy cp G [-7r,7r]^. This 
assumption is satisfied in all the test cases in Section V. 
Then, we are interested in solving the following problem: 
given B, /3, 

|2 



mm 



subject to 



where || • || denotes the Euclidean norm. This is the stan- 
dard linear least-squares estimation where /3 represents 
an observed signal that is to be projected onto the range 
space of B in order to minimize the normed error 
It admits a unique optimal solution: 



= [B^B) ^B^I3 

I-B{B^By^B^ 



(15) 
(16) 



B. Arbitrary network of phase shifters 

More generally, consider a network with phase shifters 
on an arbitrary subset of links. Given a relaxed solution 
y, under what condition does there exists a 6 such that 
the inverse projection ho{y) is a branch flow solution 
in X? If there is a spanning tree T such that all links 
outside T have phase shifters, then Theorem 1 says that 
such a 9 always exists, with an appropriate choice of 
phase shifter angles on non-tree links. Suppose no such 
spanning tree exists, i.e., given any spanning tree T, there 
is a set E_\_f Q E \ Et of links that contain no phase 
shifters. Let B^> and denote the submatrix of B 
and sub vector of /3, respective, corresponding to these 
links. Then a necessary and sufficient condition for angle 
recovery is: there exists a spanning tree T such that the 
associated B^> and satisfy: 



This condition reduces to (8) if there are no phase shifters 
in the network {E^' ^ E\ Et) and is always satisfied if 
every link outside any spanning tree has a phase shifter 
{E^f = 0). It requires that the phase angle differences 
implied by y sum to zero around any loop that contains 
no phase shifter (c.f. Theorem 2(1) and Remark 4 of 
[2]). After such a T is identified, the above two methods 
can be used to compute the required phase shifts. 

C. Other properties 

We close this section by discussing two properties of 

</.. 

First, the voltage angles are 9 = Bj}{I5t — 0r) pro- 
jected onto [— TT, tt]^, and the angle recovery condition 
(8) becomes 



(18) 



which can always be satisfied by appropriate (nonunique) 
choices of (p. A similar argument to the proof of Theorem 
2(2) leads to the following interpretation of (18). For any 
link (i, j) G E, (14) says that the phase angle difference 
from node i to node j is /3ij and consists of the voltage 
angle difference 9i — 9j = /3ij — (pij and the phase shifter 
angle (pij. Fix any link (i, j) G E\Et not in tree T. The 
left-hand side [fi^fi^^(/3r — 0r)]^^ of (18) represents 
the sum of the voltage angle differences from node i to 
node j along the unique path in T, not including the 
phase shifter angles along the path. This must be equal 
to the voltage angle difference — across (the 
non-tree) link (i, j), not including the phase shifter angle 
across (i, j). Therefore (18) has the same interpretation 
as before that the voltage angle differences sum to zero 
around any cycle, though, with phase shifters, the voltage 
angle differences are now /3^j — (pij instead of Pij. This 
in particular leads to a relationship between any two 
solutions (0*,0*) and (^,0) of (14). 

In particular, let (0*, 0*) be the solution in Theorem 1 
where 0* ^ T-^, and (^,0) any other solution. Then 
applying (18) to both 0* and leads to a relation 
between them on every basis cycle. Specifically, let 
i ^ j be a link not in the spanning tree T, let T(0 ^ k) 
be the unique path in T from node to any node k. Then 
for each link i j in E that is not in T, we have 



[0* 



(17) 



{k,i)eT{(y^j) {k,i)eT{o^i) 
{k,i)eT{(y^j) {k,i)eT{o^i) 



Second, Theorem 1 implies that given any relaxed 
solution there exists a G T-^ such that its inverse 
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projection x := he{y) is a branch flow solution, i.e., 
satisfies (9)-(ll). We now give an alternative 
direct construction of such a solution {x^ (j)) from any 
given branch flow solution x and phase shifter setting 
that may have nonzero angles on some links in T. 
It exhibits how the effect of phase shifters in a tree is 
equivalent to changes in voltage angles. 

Fix any spanning tree T. Given any partition 
0^ = [^T 0±] with respect to T. Define a G [— tt, tt]^ by 



Bt oi — 



or a := Brj 



bx. Then define the mapping 



by 



Vi : — V^e'^', /jj : — lije'^', Sij : — Sij (19) 



and 



if G 



(ai-aj) if {i,j)eE\ET 



(20) 



i.e., (/) is nonzero only on non-tree links. It can be verified 
that ai - aj = J2eeT(i^j) 4>e where T{i ^ j) is the 
unique path in tree T from node i to node j. Note that 
\Vi\ = |t^|, = \Iij\ and S = S. Hence if h{x) is a 
relaxed branch flow solution, so is h{x). Moreover, the 
effect of phase shifters in T is equivalent to adding ai 
to the phases of Vi and lij. 

Theorem 2: Fix any tree T. If (x, 0) is a solution of 
(9)-(ll), so is (x, 0) = ^(x, 0) defined in (19)-(20). 

Proo/- 1 Since |V-| = \Iij\ = |4-| and 5 = 5, 
(x,0) satisfies (lO)-(ll). For any link G E't in 

tree T, (19)-(20) imply 



= (v. 



where the second equality follows from Bt ol — ^t- For 
any link (i, j) ^E\Et not in T, (19)-(20) imply 



- Vj e-'"^'^ = (Vi - Vj e-'^'^^ e^^^ 



But (Vi — Vje = lij since (x,(^) satisfies (9). 

Therefore Vi — Vj e~^^'^ = lij, i.e., (x,(/)) satisfies (9) 
on every link. ■ 

IV. OFF IN CONVEXIFIED NETWORK 

Theorem 1 suggests using phase shifters to convexify 
a mesh network so that any solution of OFF-ar can 
be mapped into an optimal solution of OFF of the 
convexified network. Convexification thus modifies a 

less direct proof is to observe that (18) and a = B^^^t imply 
B^B^^^Pt = f3±- (0_L - B^a) = I3±-(I)± 
which means (x, 0) satisfies (14). 



NP-hard problem into a simple problem without loss in 
optimality; moreover this requires an one-time deploy- 
ment cost for subsequent operational simplicity, as we 
now show. 

We will compare four OFF problems: the original OFF 
defined in [2]: 
OPF: 

min / ( /i(x), s ) 

x,s \ J 

subject to X G X, {S^ sq, G S 

the relaxed OFF-ar defined in [2]: 
OFF-ar: 

min / ( /z(x), s ] 

x,s \ J 

subject to X G Y, [S^ sq, G S 

the following problem where there is a phase shifter on 

every line (0 G [— 7r,7r]^): 

OPF-ps: 

f{h{x),s) 
subject to X G X, {S^v^ sq, s) e S> 



mm 



and the problem where, given any spanning tree T, there 

are phase shifters only outside T: 

OPF-ps(T): 



min / ( /i(x), s 

X,S,(f) \ 

subject to X G Xt, {S, v, sq, G S, G 

Let the optimal values of OFF, OFF-ar, OFF-ps, and 
OFF-ps(T) be /*, far, fps, and /r respectively. 

Theorem 1 implies that X C Y = X = X^^ for any 
spanning tree T. Hence we have 

Corollary 3: For any spanning tree T, > far = 
fps = /r, with equality if there is a solution {yar^Sar) 
of OFF-ar that satisfies (8). 

Corollary 3 has several important implications: 

1) Theorem 1 in [2] implies that we can solve OFF- 
ar efficiently through conic relaxation to obtain a 
relaxed solution (yar^Sar)- An optimal solution of 
OFF may or may not be recoverable from it. If 
yar satisfies the angle recovery condition (8) with 
respect to Sar, then Theorem 2 in [2] guarantees a 
unique 6^ G [— tt, tt]^ such that the inverse projec- 
tion {ho^{yar)^Sar) is indeed optimal for OPF. 

2) In this case. Corollary 3 implies that adding any 
phase shifters to the network cannot further reduce 
the cost since /* = far = fps- 

3) If (8) is not satisfied, then yar ^ h{X) and there 
is no inverse projection that can recover an opti- 
mal solution of OFF from (yar^^ar)- In this case. 
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/* > far- Theorem 1 implies that if we allow 
phase shifters, we can always attain far = fps 
with the relaxed solution (^^r? Sar), with potentially 
strict improvement over the network without phase 
shifters (when > far)- 
4) Moreover, Corollary 3 implies that such benefit can 
be achieved with phase shifters only in branches 
outside an arbitrary spanning tree T. 

V. Simulations 

For radial networks, results in Part I (Theorem 4) 
guarantees that both the angle relaxation and the conic 
relaxation are exact. For mesh networks, the angle relax- 
ation may be inexact and phase shifters may be needed 
to implement a solution of the conic relaxation. We now 
explore through numerical experiments the following 
questions: 

• How many phase shifters are typically needed to 
convexify a mesh network? 

• What are typical phase shifter angles to implement 
an optimal solution for the convexified network? 

Test cases. We explore these questions using the IEEE 
benchmark systems with 14, 30, 57, 118 and 300 buses, 
as well as a 39-bus model of a New England power 
system and two models of a Polish power system with 
2,383 and 2,737 buses. The data for all the test cases 
were extracted from the library of built in models of 
the MATPOWER toolbox [5] in Matlab. The test cases 
involve constraints on bus voltages as well as limits on 
the real and reactive power generation at every generator 
bus. The New England and the Polish power systems also 
involve MVA limits on branch power flows. All these 
systems are mesh networks, but very sparse. 

Objectives. We solve the test cases for two scenarios: 

• Loss minimization. In this scenario, the objective is 
to minimize the total active power loss of the circuit 
given constant load values, which is equivalent to 
minimizing the total active power generation. The 
results are shown in Table I. 

• Loadability maximization. In this scenario, the ob- 
jective is to determine the maximum possible load 
increase in the system while satisfying the genera- 
tion, voltage and line constraints. We have assumed 
all real and reactive loads grow uniformly, i.e., 
by a constant multiplicative factor called the max 
loadability in Table II. 

Solution methods. We use the "SEDUMI" solver in 
Matlab [6]. We first solved the SOCP relaxation OPF-cr 
for a solution (^, s) of OPF-ar. In all test cases, equality 
was attained at optimality for the second-order cone 
constraint, and hence OPF-cr was exact, as Theorem 1 in 



[2] would suggest. Recall however that the load values 
were constants in all the test cases. Even though this 
violated our condition that there are no upper bounds 
on the loads OPF-cr turned out to be exact with respect 
to OPF-ar in all cases. This confirms that the no-upper- 
bound condition is sufficient but not necessary for the 
conic relaxation to be exact. 

Using the solution s) of OPF-ar, we checked if 
the angle recovery condition (8) was satisfied. In all test 
cases, the angle recovery condition failed and hence no 
{ho{y)^s) was feasible for OPF without phase shifters. 
We computed the phase shifter angles (p using both meth- 
ods explained in Section III-A and the corresponding 
unique {he{y)^s) that was an optimal solution of OPF 
for the convexified network. For the first method that 
minimizes the number of required phase shifters, to place 
the phase shifters, we have used a minimum spanning 
tree of the network where the weights on the lines are 
their reactance values. 

In Tables I and II, we report the number m—n of phase 
shifters potentially required, the number of active phase 
shifters (i.e., those with a phase angles greater than 0.1°), 
and the range of the phase angles at optimality using both 
methods. In Table II, we also report the simulation time 
on an Intel 1.8 GHz Core 15 CPU. 

We report the optimal objective values of OPF with 
and without phase shifters in Tables I and II. The optimal 
values of OPF without phase shifters were obtained 
by implementing the SDP formulation and relaxation 
proposed in [7] for solving OPF. In all test cases, the 
solution matrix was of rank one and hence the SDP 
relaxation was exact. Therefore the values reported here 
are indeed optimal for OPF. 

The SDP relaxation requires the addition of small 
resistances (10~^ pu) to every link that has a zero 
resistance in the original model, as suggested in [8]. This 
addition is, on the other hand, not required for the SOCP 
relaxation: OPF-cr is tight with respect to OPF-ar with 
or without this addition. For comparison, we report the 
results where the same resistances are added for both the 
SDP and SOCP relaxations. 

Summary. From Tables I and II: 

1) Across all test cases, the convexified networks 
have higher performance (lower minimum loss and 
higher maximum loadability) than the original net- 
works. More important than the modest perfor- 
mance improvement, convexification is design for 
simplicity: it guarantees an efficient solution for 
optimal power flow. 

2) The networks are (mesh but) very sparse, with the 
ratios m/(n + 1) of the number of lines to the 
number of buses varying from 1.2 to 1.6 (Table I). 
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No PS 


With phase shifters (PS) 


Test cases 


# links 


Min loss 


Min loss 


# required PS 


# active PS 


Min #PS (°) 


Min n 




(m) 


(OPE, MW) 


(OPF-cr, MW) 


(m 


-n) 


>0.1° 


[0min 7 0max] 


[0min 7 0max] 


IEEE 14-Bus 


20 


0.546 


0.545 


1 


(35%) 


2 (10%) 


[-2.09, 0.58] 


[ -0.63, 0.12] 


IEEE 30-Bus 


41 


1.372 


1.239 


12 


(29%) 


3 (7%) 


[-0.20, 4.47] 


[-0.95, 0.65] 


IEEE 57-Bus 


80 


11.302 


10.910 


24 


(30%) 


19 (24%) 


[-3.47, 3.15] 


[-0.99, 0.99] 


IEEE 118-Bus 


186 


9.232 


8.728 


69 


(37%) 


36 (19%) 


[-1.95, 2.03] 


[-0.81, 0.31] 


IEEE 300-Bus 


411 


211.871 


197.387 


112 


(27%) 


101 (25%) 


[-13.3, 9.40] 


[-3.96, 2.85 ] 


New England 39-Bus 


46 


29.915 


28.901 


8 


(17%) 


7 (15%) 


[-0.26, 1.83] 


[-0.33, 0.33] 


Polish (case2383wp) 


2,896 


433.019 


385.894 


514 


(18%) 


373 (13%) 


[-19.9, 16.8] 


[-3.07, 3.23] 


Polish (case2737sop) 


3,506 


130.145 


109.905 


770 


(22%) 


395 (11% ) 


[-10.9, 11.9] 


[ -1.23, 2.36] 



TABLE I: Loss minimization. Min loss without phase shifters (PS) was computed using SDP relaxation of OFF; 
min loss with phase shifters was computed using SOCF relaxations OFF-cr of OFF-ar. The "(%)" indicates the 
number of FS as a percentage of #links. 





No PS 


With phase shifters (PS) 


Test cases 


Max load. 


Max load. 


# required PS 


# active PS 


Min #PS (°) 


Min ||(/) 


r n 


Simu. time 




(OPF) 


(OPF-cr) 


(m 


-n) 


\M 


> 0.1° 


[0min 7 0max] 


[0min7 


^max] 


(seconds) 


IEEE 14-Bus 


195.0% 


195.2% 


1 


(35%) 


6 


(30%) 


[-0.51, 1.35] 


[-0.28, 


0.23] 


1.92 


IEEE 30-Bus 


156.7% 


158.7% 


12 


(29%) 


9 


(22%) 


[-0.42, 12.4] 


[-2.68, 


1.86] 


3.86 


IEEE 57-Bus 


108.2% 


118.3% 


24 


(30%) 


24 


(30%) 


[-13.1, 23.2] 


[-4.12, 


4.12] 


7.13 


IEEE 118-Bus 


203.7% 


204.9% 


69 


(37%) 


64 


(34%) 


[-8.47, 17.6] 


[-4.61, 


5.36] 


15.96 


IEEE 300-Bus 


106.8% 


112.8% 


112 


(27%) 


103 


(25%) 


[-15.0, 16.5] 


[-4.28, 


6.31] 


34.6 


New England 39-Bus 


109.1% 


117.0% 


8 


(17%) 


5 


(11%) 


[-1.02, 1.28] 


[-0.28, 


0.18] 


2.82 


Polish (case2383wp) 


101.4% 


106.6% 


514 


(18%) 


435 


(15%) 


[-19.6, 19.4] 


[-4.06, 


4.32] 


434.5 


Polish (case2737sop) 


127.6% 


132.5% 


770 


(22%) 


420 


(12%) 


[-13.9, 17.1] 


[-2.07, 


3.62] 


483.7 



TABLE II: Loadability maximization. IVIax loadability without phase shifters (PS) was computed using SDP 
relaxation of OPF; max loadability with phase shifters was computed using SOCP relaxations OPF-cr of OPF- 
ar. The "(%)" indicates the number of PS as a percentage of #links. 



The numbers m — n of phase shifters potentially 
required on every link outside a spanning tree for 
convexification vary from 17% of the numbers m 
of links to 37%. 

3) The numbers of active phase shifters in the test 
cases vary from 7% of the numbers m of links to 
25% for loss minimization, and 11% to 34% for 
loadability maximization. The phase angles required 
at optimality is no more than 20° in magnitude with 
the minimum number of phase shifters. With the 
maximum number of phase shifters, the range of 
the phase angles is much smaller (no more than 
5°). 

4) The simulation times range from a few sec to 
mins. This is much faster than SDP relaxation. 
Furthermore the simulation time appears linear with 
the network size. 

VI. Conclusion 

We have presented a branch flow model and demon- 
strated how it can be used for the analysis and opti- 
mization of mesh as well as radial networks. Our results 
confirm that radial networks are computationally much 
simpler than mesh networks and we should exploit this 



advantage whenever we can. For mesh networks, we 
have proposed a simple way to convexify them using 
phase shifters that will render them computationally 
as simple as radial networks for power flow solution 
and optimization. In other words, the addition of phase 
shifters convert a nonconvex problem into a different, 
simpler problem. 

We have proposed a solution strategy for OPF that 
consists of two steps: 

1) Compute a relaxed solution of OPF-ar by solving 
its conic relaxation OPF-cr. 

2) Recover from a relaxed solution an optimal solu- 
tion of the original OPF using an angle recovery 
algorithm. 

We have proved that, for radial networks, both steps 
are always exact, provided there are no upper bounds 
on loads, so this strategy guarantees a globally optimal 
solution. For mesh networks the angle recovery condition 
may not hold but can be used to check if a given relaxed 
solution is globally optimal. 

Since power networks in practice are very sparse, 
the number of required phase shifters may be relatively 
small. IVIoreover, the placement of these phase shifters 
depends only on network topology, but not on power 



flows, generations, loads, or operating constraints. There- 
fore only one-time deployment cost is required to achieve 
the subsequent simplicity in network operation, market 
operation, and investment decisions. 
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